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(N : ABSTRACT 

m 

^ I We obtain an unbiased estimate of the global Hubble constant Hq in the volume of 

cz < 12000 km s-i in the region of Pisces-Perseus. The Tully-Fisher (TF) relation is 
CSl ' applied to a magnitude-limited sample of 441 spirals selected from the Arecibo 21 cm 

^ , catalog. We improve the photometry data of the previous TF study by Ichikawa and 

I Fukugita (1992) by using our original surface photometry data and local calibrators 

' which are newly made available in the literature. The photometry data were calibrated 

• with CCD observations and we achieve 0.13 mag for a photometric internal error. 

We use a maximum likelihood method for the TF analysis. Monte-Carlo simulations 



O I demonstrate that our method reproduces a given Hq at the 95% confidence level. By 

' applying the method to our sample galaxies, we obtain the unbiased global Hubble 

Ci ■ constant Hq = 65 it 2 kms~^Mpc~^ ; the first and the second terms represent the 

internal random error and the external uppermost and lowermost systematic errors, 
^ I respectively. We also find a good agreement for our Hq with those recently obtained 

^ ' via Cepheid observation in the local Universe and the TF relation and supernovae 

applied to a spatial volume comparable with or larger than ours. Hubble velocities of 
the spirals inferred from our Hq show no significant systematic difference from those 
given in the Mark III catalog. The same analysis for Hq is carried out using r-band 
photometry data of the Pisces-Perseus region given by Willick et al. (1997). We obtain 
a global Hq which is consistent with that obtained from the i?-band analysis. A bulk 
motion in the Pisces-Perseus region is briefly discussed, based on our calibration of 
Hq. The -B-band intrinsic TF scatter is too large to allow any determination of bulk 
motion. However, our r-band TF analysis supports the notion of a coherent streaming 
motion of the Pisces-Perseus ridge with ~ — 200kms~^ with respect to the CMB, in 
agreement with most modern studies. 



Subject headings: distance scale — galaxies: clusters : individual (Pisces-Perseus) — 
galaxies: distances and redshifts 
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1. Introduction 

The global value of the Hubble constant Hq is determined by investigating a Hubble flow of 
a large number of galaxies distributed in a volume larger than the scale of inhomogeneous galaxy 
distribution in the Universe (~ 100/i~^Mpc, h = i7o/100kms~^Mpc~^). Ichikawa & Fukugita 
(1992, hereafter IF92) determined the global Hq by applying the Tully-Fisher (TF) relation (Tully 
&; Fisher 1977) to the Arecibo HI 21cm sample ( Giovanelli and Haynes 1985| , piovanelli et al.| 



1986, Giovanelli and Haynes 1989| ). The sample extends in the area of 950 square-degrees in the 



region of Pisces-Perseus (22 < 01959 < 4 , -|-22° < ^igso < +33°) and includes all the field and 
cluster galaxies complete down to = 15.5 mag originated from CGCG (Zwicky et al. 1961-68). 
Because of the large area and the deep limiting magnitude, the sample IF92 used contains about 
400 spirals at cz < 12000 km s~^. This large sample enabled them to examine an effect of biases in 
the TF analysis. 

In this paper, we re-determine Hq using the Arecibo HI sample. The differences from the 
study of IF92 are as follows. First, the imaging data were extracted from our photographic plate 
material via modern image reduction techniques. IF92 used magnitudes and axial ratios which 
were transformed into a homogeneous system of RC2 (de Vaucouleurs, de Vaucouleurs & Corwin 
1976) from CGCG, UGC (Nilson 1973), and the Arecibo catalog. Second, we use a new set of 
local calibrators for TF calibration. Their distances are all calibrated with Cepheid observations. 
Third, our TF analysis uses a technique of maximum likelihood method. Fourth, our TF analysis 
is also applied to r-band CCD data of "w91pp" in the Mark III catalog (Wilhck et al. 1997). The 
sample was originated from UGC and hence complete down to = 14.5 mag, shallower than 
ours by 1 mag. Nevertheless, the data allow us to carry out a useful comparison with our result 
because the r-band TF relation is tighter than that of S-band. Finally, we discuss, based on our 
result of the global Hq, the bulk motion of galaxies in the region which has been studies by many 
authors (WiUick 1990, Han and Mould 1992, Courteau et al. 1993, Hudson et al. 1997). 

Our adoption of photographic data may deserve a comment. Motivation of the present study 
is to determine Hq in the Pisces-Perseus region where a bulk motion of galaxies has been attracted 
much attention (Wilhck 1991 ; Han and Mould 1992 ; Courteau et al. 1993 ; Hudson et al. 1997). 
The bulk motion has been usually computed using relative distances of galaxies with respect to a 
reference point of the Hubble flow. In this procedure, the motion is derived independently from 
Hq. In this paper, on the other hand, we attempt to determine the global Hq and then discuss the 
bulk flow using the absolute value of Hq. 

To achieve the goal of our study, our galaxy sample should fulfill three conditions. First, it 
should be large in terms of the spatial volume so that the global Hq is determined unaffected 
by the bulk flow. The volume is required to be larger than the coherent scale of the flow 
{cz ~ 6000 km s~^). Second, our sample should be also large in terms of the number of galaxies 
and should not be sparse so that we can evaluate the effect of a bias correctly. Third, it is of course 
highly desired that photometry data are homogeneous so as to eliminate systematic errors which 
creep in a data set taken from different compilations. Considering these conditions, we carried 
out photographic observations and obtained homogeneous photometry data for the Pisces-Perseus 
region. The galaxy sample is a magnitude limited one containing all the CGCG/UGC galaxies, and 
is distributed up to cz ~ 12000 km s~^. The sample is therefore larger than a Pisces- Perseus fleld 
sample "w91pp" in the Mark III catalog (Willick et al. 1997). A recent study of Hq by Giovanelli 
et al. (1997a) covered the entire sky using 24 clusters. Although the Pisces-Perseus region was 
included in their research, their sample in the region is limited to five clusters including about 150 
galaxies. Further, the accuracy of our photographic photometry is as good as 0.13 mag ( |Watana"be 



1996). It is worth stressing that careful photographic photometry do achieve 0.1-0.2 mag accuracy 
(e.g., Kodaira et al. 1990 ; Yasuda, Okamura and Fukugita 1995). Thus, our photographic 
observations offer the photometry data which are most suitable for our study. 

Outline of this paper is as follows. The data used for the present analysis are described in 
§^. The maximum-likelihood method of the TF analysis is presented in Calibration of the TF 
relation is given in §|^. Reliability of the present method is examined by Monte-Carlo simulations 
in §|5[ The method is applied to our i?-band sample and the unbiased estimate of the global 
Ho is determined in §^. A comparative analysis is carried out using "w91pp" data of Mark III 
compilation in §0. We discuss the present result in §|8[ 



2. Data 

2.1. Observations and Data Reduction 

Observations and data reduction are briefly described here. Details will be given in a 
forthcoming paper. 

Photographic observations of the region of Pisces-Perseus {22^ < 01950 ^ and 
+22° < (5i95o < -|-33°) were carried out in B (Kodak Ha-O emulsion plus Schott GG385 glass 
filter) using the 105 cm Schmidt telescope at the Kiso Observatory, the University of Tokyo. Most 
of the observations were conducted in 1991 and in 1992 and we obtained thirty-eight plates. 
We measured and reduced all the 1524 CGCG and UGC galaxies in the region. These 1524 
galaxies are referred to as the photometry sample. The galaxies were scanned with a PDS2020 
microdensitometer at the Kiso Observatory. We then reduced the images in a conventional manner 
for Schmidt plate data (e.g., Kodaira et al. 1990) using the software packages IRAfI and SPIRAL 
(Ichikawa et al. 1987; Hamabe and Ichikawa 1991). CCD observations were also carried out in B 
and in V for seventy spirals selected from the photometry sample to calibrate the photographic 
data. The CCD data were calibrated with photometric standard stars (Landolt 1992) of various 
colors observed at various zenith distances to correct atmospheric extinction and to transform our 
photographic i?-band system to the Johnson's. For seven plates where no CCD observation was 
available, the flux calibration was made with photoelectric aperture photometry data of Longo 
and de Vaucouleurs (1983). The flux calibration of the photographic data using the CCD and the 
photoelectric aperture photometry data followed the conventional method of Kodaira et al. (1990). 
Aperture magnitudes were measured with a series of concentric circular apertures throughout 
the present study. Total magnitudes Bt were obtained by a sliding fit of aperture-magnitude 
curve (hereafter a growth curve) to their templates given in RC3 (de Vaucouleurs et al. 1991). 
Isophotal magnitudes -625 and major-to-minor axial ratios R25 were measured at the isophote of 
25 mag arcsec"^ in the surface brightness. Internal errors were 0.13 mag for Bt and 0.06 dex for 
logi?25 ([Watanabe 199^ ). 



Spiral galaxies are selected from the photometry sample and used for the present TF analysis. 
The procedure of the selection is described in §127 



'^IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science 
Foundation. 
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2.2. Comparisons with Other Data 

We examine external consistency by comparing our data with those in RC3 (N=79) and 
Cornell et al. (1987) (N=8). All overlapping systems were rejected from the comparisons. We find 
possible errors in Bt of RC3 for PGC5473 (CGCG 502-87, UGC 1045) and in Bt of RC3 and in 
B25 of Cornell et al. (1987) for PGC3133 (Watanabe 1996). These galaxies are further rejected. 
The comparisons are presented in Fig.|^ and Fig.^. From a regression analysis for RC3 data, we 
obtain BT(Kiso) = 5r(RC3) + 0.02 mag {a = 0.16 mag) and log i?25(Kiso) = logi?25(RC3) + 0.00 
(fj = 0.07 dex). If we assume that the half of u = 0.16 mag is ascribed to the error in RC3, 
our Bt is supposed to be as accurate as 0.11 mag. This is compatible with our internal error 
of 0.13 mag for Bt- Rms errors of Bt between RC3 data and ours are 0.17 mag (N=46) for 
galaxies calibrated with our CCD observation while 0.14 mag (N=32) for those calibrated with 
the aperture photometry data by Longo and de Vaucouleurs (1983). The difference between 
these two values is insignificant, and we do not think the calibration by the aperture photometry 
data is poorer than that by the CCD data. A regression analysis for Cornell's data leads to 
S25(Kiso) = B25(Cornen) + 0.01 mag {a = 0.01 mag) and log i?25(Kiso) = log i225(Cornen) - 0.01 
(fj = 0.04 dex). These good agreements with their CCD data demonstrate reliability of our 
photographic photometry data. 



2.3. HI Data 



HI 21 cm data are taken from the Arecibo HI catalog ( Giovanelli and Haynes 1985 , Giovanelli 



et al. 1986 , piovanelli and Haynes 1989 ), where line-profile widths are defined at the 50% level of 

50 



the line mean intensity. The line widths W^q'^ and heliocentric recession velocities Vh both have 



an accuracy of about lOkms . 

In earlier TF calibration, the 21 cm line-profile width has been defined with two schemes. One 
refers to the width W^q"^ as used in the Arecibo catalog, while the other adopted a width ^20^ 
measured at the 20% level of the line peak intensity. We adopt the latter definition ^^20^ because 
most of the earlier TF calibrations have been made with this scheme (e.g., Aaronson et al. 1986, 
Bottinelh et al. 1987, Pierce and Tully 1988, Wilhck 1990, Pierce & Tully 1992, IF92). Adopting 
this schemes allows us to compare our TF calibration easily with those of the previous studies. 
We convert W§^' to W^^' using an empirical formula given by Willick (1991) ; 

log VF2°o' - 2.5 = 0.933 (log W^l^' - 2.5) + 0.019 . (1) 

The relation has a dispersion of 0.02 dex. The uncertainty caused by this dispersion is less than 
0.1 mag which is negligible compared with the intrinsic dispersion of the TF relation, > 0.3 mag. 

We determine the global Hq with reference to a rest frame of the cosmic microwave background 
(CMB) radiation. The heliocentric velocity V/j is accordingly transformed to Vc measured in the 
frame using CMB dipole anisotropy attributed to the solar motion of f = 365kms~^ toward 
{l,b) = (265°, +48°) relative to the CMB radiation field (Smoot et al. 1991). Let Ig and bg be 
the galactic longitude and latitude of a particular galaxy. Then, the transformation is expressed 
as Vc = Vh + Vx + Vy + Vz where Vx = v cos b cos / cos bg cos Ig, Vy = v cos b sin / cos bg sin Ig, and 
Vz = V s'mbs'mbg, respectively. 



The Arecibo catalog mainly refers to UGC for galaxy morphological types. For galaxies for 



- 5 - 



which the type was not available in UGC, authors of the Arecibo catalog determined the type 
from visual inspection of galaxy images on Palomar Schmidt plates or on Palomar Sky Survey 
prints. We take their morphological type with a notation Tgh- 



2.4. Correction to the Observational Parameters 

2.4-1- Apparent magnitudes 

The observed magnitude Bj- is corrected for Galactic and internal extinction and K-dimming. 
The Galactic extinction Aq is taken from a machine-readable extinction map made by David 
Burstein (hereafter referred to as BH). The internal extinction Ai and the K-dimming K are 
both evaluated following a prescription given in RC3. Corrected magnitude Bj^ is then given by 
B^ = BT-AG-Ai-K. 

There are various schemes used for evaluation of Ac (e.g., BH, RC2, RSA [Sandage and 
Tammann 1981]) and Ai (e.g., RC3, RC2, RSA, Bothun et al. [1985]). Systematic offset among 
these different schemes will be taken into account as a source of an error for Hq (§|6.1D. 



2-4-2- Line-profile widths 
The line-profile width VF^o'^ corrected for a projection and a cosmological broadening effects 

with 

W,o = ,^ 7? ■ . , (2) 

(1 + z) Sim 

where inclination i of a spiral is estimated from i?25 via a formula 



(3) 



Dependence of the intrinsic major-to-minor axial ratio Rq on a galaxy type was investigated 
by some authors (e.g., Bottinelli et al. 1983, Fouque et al. 1990, Giovanelli et al. 1997b| ). Some 




authors also favors a low value Rq < 0.2 (Giovanelli et al. 1994 ; Giovanelli et al. 1997b ; Willick et 
al. 1997). However, we use a conservative value Rq^ = 0.2 for all spirals following most of earlier 
analysis (Pierce and Tully 1988 ; Fukugita et al. 1991 ; Freudling, Martel and Haynes 1991 ; IF92 
; Han 1992 ; Pierce and Tully 1992 ; Bernstein et al. 1994 ; Bureau, Mould and Staveley-Smith 
1996). We find that adopting a low value of Rq^ = 0.13 or introducing its type dependence results 
in an insignificant change in Hq by ~ Ikms^^Mpc^^ in our analysis. 

Giovanelli et al. (1997b) proposed a different formula for the correction of the line-profile 
width W^Q. A similar formula could be worked out for W2o- However, we did not pursue this way. 
The ratio of their observed width W to the corrected width Wi is almost constant ; the mean 
of W/Wi is 1.07 with a scatter of only 0.07. It is reasonable to assume that the ratio would be 
nearly constant in case such a revision is made for W2o- Further, such a revision would change 
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W20 of local calibrators as well as sample galaxies. Then, the effect of such a revision is likely to 
be canceled out. 



2.5. Sample Selection 

We select sample galaxies for our TF analysis according to selection criteria as follows; 1) 
Morphological type Sa-Sd, i.e., 3 < Tgh < 8 and Tgh = 12. This criterion is adopted to exclude 
SO and very late-type spirals unsuitable for the TF analysis. Galaxies classified to Tqh = 12 
are all regarded as Sb in the present analysis. 2) Inclination i > 45°. For galaxies with i < 45°, 
an uncertainty of the TF magnitude caused by an error of Alogi?25 = 0.06 (Watanabe 1996) 



combined with a slope of 6.14 for our TF calibration (§4.1) exceeds 0.4 mag, which is no more 
negligible compared with the intrinsic dispersion of the TF relation. 3) Recession velocity 
Vc > 1000 km s~^. This criterion is adopted to avoid disturbance by a local peculiar velocity field 
with amplitude of several hundreds kms~^. 4) Line-profile width 2.30 < logW20 < 2.75. As 



discussed in §4.1, the lower and the upper boundaries are originated from the range of log 1^20 for 



the local calibrators we adopted. 

These criteria reduced the photometry sample to 449 galaxies (Table ||). These are hereafter 
referred to as the spiral sample. The TF analysis is carried out using 441 galaxies among them 
(§^. The eight galaxies are discarded from the analysis because of their large deviation in the 



TF magnitude; this is discussed in §0. A wedge diagram for the spiral sample is shown in Fig.^, 
which illustrates that the sample covers a spatial volume up to ~ 12000 km s^^. Number 
distribution as a function of Bt is shown in Fig.^. The number ratio of the spiral sample to the 
photometry sample is almost constant over 13 < Bt < 17 mag. The spiral sample is therefore 
supposed to be a fair subset of the photometry sample. 



3. A MAXIMUM LIKELIHOOD TECHNIQUE 

If we obtain a TF distance rj-p using Bj, and a TF magnitude Mtf as vtf = 10^^t~~^^tf+5)/5 
the most straightforward estimate of Hq may be a ratio Vc/ttf- We hereafter call it as the Hubble 
ratio. For our spiral sample, however, a sample average of the Hubble ratios is not a correct 
estimate of the true global Hq. It is subject to a systematic error because the sample galaxies are 
biased toward intrinsically bright ones. This bias is originated from the facts that the intrinsic 
dispersion cjtf of the TF relation is non-zero and that the spiral sample is a magnitude-limited 
sample. Consequently, a systematic error occurs in the sample average of the Hubble ratio. In this 
paper, we refer to this bias as a sampling incompleteness bias. We also take the Malmquist bias 
(Malmquist 1922) into account in the present analysis. The Malmquist bias is caused by a non-zero 
intrinsic dispersion of the TF relation combined with a conical shape of an observed volume for a 
unit solid angle. A detailed review for the two biases is found in Strauss and Willick (1995). Note 
that an error caused by the Malmquist bias is independent of a limiting magnitude, while an error 
caused by the sampling incompleteness bias depends basically on the limiting magnitude. 

To determine a global Hq which is not affected by such systematic errors, we use a maximum 
likelihood method in the TF analysis. A similar technique was used by Han and Mould (1990) 
for an analysis of a Local Supercluster velocity field. The likelihood function is computed using a 
probability distribution of W20 rather than Bt since it is less affected by the sample incompleteness 
(e.g., Schechter 1980; Freudling et al. 1995, and references therein). This is equivalent to the 
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"inverse method" discussed by Strauss and Willick (1995). 

3.1. Formulation 

Let us consider a probability P*-*-* that the i-th spiral galaxy has a line-profile width 1^20"*''*^ 
under the condition that it has a true absolute magnitude M^*\ an inclination and a redshift 
z^*^. Note that M^^ is an unknown parameter. An intrinsic line-profile width is obtained 
from W^Q^^*"*, and z^^^ with equations and (|3|). A TF magnitude corresponding to 

is then obtained by the TF relation, 

M^l = A + B (log W^S - 2-5) . (4) 

(i) 

On the other hand, the Bubble's law gives an estimate for , i.e., 

m};) = - 5 log VP + 5 log /i - 15 - (5 , (5) 

where is given in unit of kms~^. The additional term 5 is a correction for a systematic error 
caused by the Malmquist bias. If galaxies are uniformly distributed in space, the systematic error 
is (5 = (3 In 10/5) = 1.38(j^ , where cr is a random distance error expressed in term of a distance 
modulus. In the present analysis, the measure of the distance is Vc (equation (^)). Accordingly, 
the random distance error comes from a random error AVc, which is evaluated as a square root of 
a quadratic sum of a line-of-sight component of a random peculiar velocity of galaxies Crand 

and 

an observational error of the recession velocity AV/i, 



(6) 



Using AV^, we obtain a = 51oge • {/S.Vc/Vc'^^). The systematic error 5 caused by the homogeneous 
Malmquist bias is therefore given by 



(5 = 15 log e 



AK 



(7) 



An effect of a spatial inhomogeneity of galaxy distribution is discussed later in §5.2 



We assume that an intrinsic broadening of the TF relation follows a Gaussian function with a 

and 



dispersion utf- Then, using mI^, and we obtain the probability 7^^*^ in a following form: 



1 



exp 



{mI^p - M 



«^2 
B ) 



(8) 



where cT(j) includes cttf and uncertainties of M^^p and of M^*\ i.e.. 



r«^2■ 



(9) 
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The uncertainty /S.Mj.p is evaluated from equations and (Q) as 

where c is the velocity of light. The uncertainty AM^''', on the other hand, is computed from 
equation (^) as 

1 

2 

(11) 

With the selection criterion 2.30 < log 1^20 < 2.75 (§ ^.51) , the factor (7*^*^ is computed by 

where Mj,p, and M^p are TF magnitudes corresponding to the boundaries of log 1^20 
2.75 , respectively. 

The logarithmic likelihood function is finally obtained as 



obs(i) 



+ 



20 



.c + V, 



(0 



+ 



In 10 • A log i?25 



(i)2 
25 



(10) 



AM 



(A5, 



0W> 

T > 



+ 5 log e • 



AF, 



Vr 



(12) 
= 2.30 and 



log£ = logH^^'^ 



log e • 



(M 



(0 

TF 



M 



{i)\2 



E 



log C^^ 



(13) 



The most probable values for parameters in question are determined by maximizing log£. 

We first apply the method to local calibrators to calibrate the TF relation (jj |4.l| ). In this 
case, we let A, B in equation (^) and axp in equation (|9|) be free parameters. The method is then 
applied to our TF sample (§^) with the calibrated TF relation to determine the unbiased Hq. In 
this case, we let h in equation (^) and cjtf in equation (|9|) be free parameters. 

There are five preset parameters to specify in our method; AWJo^, CFrand, AV/i, Alogi?25, and 
ABj,. When the method is applied to local calibrators, values for AWJo^, Alogi225, and ABj. are 
taken from RC3. The parameters (Trand and AVh are unnecessary in this case because we do not 
use equation (^) but use a known distance modulus to obtain Mb for local calibrators. When the 
method is applied to our TF sample, on the other hand, following values are used in common to all 
the galaxies in the sample. AW20* ~ 20kms~^, which is estimated as a square root of a quadratic 
sum of a measurement error ~ lOkms"^ for W^q"^ (Giovanelli and Haynes 1985) and an uncertainty 
~ 15kms~^ associated with the line-profile width conversion by equation (Q). arand = 300kms~^, 
which is taken from Davis and Peebles (1983), de Lapparent, Geller and Huchra (1988) and Fisher 
et al. (1994). AVh = lOkms"^ is taken from GiovaneUi and Haynes (1985). Alogi225 = 0.06 
is taken from our error assessment (Watanabe 1996). ABj. = 0.16 mag, which is evaluated as 
a square root of a quadratic sum of our photometric uncertainty ABt = 0.13 mag (Watanabe 
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1996) and uncertainties AAq = 0.07 mag and AAi = 0.05 mag which is evahiated from Alogi?25. 
Uncertainties of these values for the five preset parameters are taken into account in an assessment 
of an external error for Hq (§|6.1|). 



4. THE TULLY-FISHER RELATION 
4.1. Calibration 

The TF relation is calibrated using local calibrators whose distances have been determined 
with Cepheid variables. These calibrators are further restricted to those which meet the same 
selection criteria as imposed on the Pisces-Perseus galaxies, i.e., a morphological type Sa-Sd and 
an inclination i > 45° (§^^). Ten spirals meet the criteria and are listed in Table Among 
the ten calibrators, NGC 224, 598, and 3031 were recently suspected to be unsuitable for the 
TF calibration because of perturbations in the internal velocity field ( piovanelli et al. 1997a| ). 
Further we find that NGC 1365 and 4639 show large difference in their inclination between RC3 
and Giovanelh et al. (1997a) {i = 59° and 48° for NGC 1365 and i = 49° and 55° for NGC 4639). 
We therefore exclude the five galaxies and use the remaining five as a fiducial set of calibrators. 
Calibration using all the ten calibrators gives only slight change in the following analysis. The 
alternative result is mentioned when necessary. 

Distance moduli /u^*-* of the calibrators are taken from the literature (Table |2|). These distances 
are all based on the LMC distance modulus of 18.5 mag (e.g., Madore and Freedman 1991). Other 
observational data are all taken from RC3. Using /u^*-* and apparent magnitudes B^^\ absolute 

(i) 

magnitudes Mj^ of the local calibrators are given by 

M« = - . (14) 

A TF diagram of the calibrators is shown in Fig.^. The range of W20 covered by the local 
calibrators, 2.30 < logW20 ^ 2.75, is already adopted as one of the selection criteria for the 
spiral sample (§ |2.5| ). Letting A, B and axF be free parameters, we apply the maximum likelihood 
method and then obtained the following calibration ; 

M^'P = -19.73 (±0.20) - 6.14 (±1.60) (log W20 - 2.5) . (15) 

The scatter cttf is 0.45 mag. If we use all the ten calibrators, we obtain A= —19.67 ± 0.15, 
B= —5.78 ± 0.95, and cttf = 0.44 mag. Our slopes are both consistent with that obtained by 
Strauss & Willick (1995) (B= —5.3) within the errors, while the values are marginally smaller 
than that given by Pierce & Tully (1992) (B= —7.48). This discrepancy is ascribed to the lack of 
our calibrators with log 1^20 < 2.30 that were used by them. 



4.2. Examination of the TF Magnitude of the Spiral Sample 

( 

In determination of Mrpp , the inclination i is the only observable derived from our 
measurements (equation [^). To examine an error in M^^^ caused by an error of our i?25, we 
plot in FigJI the residual M^j^ - Mb as a function of log i?25 ■ The absolute magnitude Mb is 
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computed from the Hubble's law (equation (^)) with Hq arbitrarily assumed to be log/i = —0.15 
[Hq = 70.8kms^^Mpc^^). Accordingly, a global offset of the plot from a horizontal line of 
Mrpp — Mb = is of little significance. We found no significant systematic trend for Mj^p — Mb 
as a function of logi?25 but found several deviant galaxies. Dispersion of the plot is 0.53 mag 
and, if we assume a Gaussian distribution for Mj,p — Mb for 449 galaxies, less than 0.5 galaxies 
fall in the range out of 3.06(t. Accordingly, deviation larger than \Mj.p — Mb\ = 1.62 mag is 
considered to be significant, and we find eight galaxies out of the boundaries. Four out of the 
eight have possible reasons for the deviation as follows. PGC 6802 has very low surface brightness 
and then unreliable Bt and -R25- PGC 10013 has a very wavy growth curve which led unreliable 
Bt- PGC 10740 is extremely faint {Bt = 17.09 mag) and hence its photometry may be unreliable. 
PGC 14033 suffers from severe contamination of HI flux of nearby galaxies ( piovanelli and Haynes 



1985). We discard these from the TF analysis. Although we found no reason for the large deviation 



to the remaining four galaxies, we further exclude them from the TF analysis and hereafter refer 
to the remaining 441 galaxies as the TF sample (Table ||). An effect of including the latter four 
galaxies is mentioned when necessary. The velocity-distance relation is shown in Fig.|7[ where 
global agreement with the Hubble's law is evident. 



5. MONTE-CARLO SIMULATIONS 

5.1. Simulated Sample 

Evaluation of the effect of biases is an important issue in this study. However, it is quite 
difficult to evaluate the effect analytically because physical properties, spatial distribution, and 
observational limits for galaxies are entangled so complicatedly. For this purpose, a Monte-Carlo 
simulation provides a straightforward and convincing test (e.g., Kolatt et al. 1996). 

We perform Monte-Carlo simulations using pseudo samples of 441 galaxies. These samples 
are made with unknown parameters given a priori as log/i = —0.15 {Hq = 70.8kms^^Mpc~^) 
and aTF = 0.40 mag. Procedure of assigning observational parameters to individual galaxies is 
described below. 

Galaxies are distributed uniformly with respect to a pure Hubble flow velocity Vh- In the 
range of 4000 <Vh < 5000kms-\ five times as many galaxies as those in the uniform background 
are added to imitate the clustering feature of the Pisces-Perseus supercluster. Galaxies in 
1000 < Vh < 4000 km s~^ are randomly discarded so that the number density becomes one third 
of that of the uniform background so as to model the low number density space in the foreground 
of the supercluster (Haynes and Giovanelli 1986). A true distance is given to the galaxies via the 
Hubble's law with Vh and log h = —0.15. We then define a velocity Vc by adding to Vh a Gaussian 
random motion with dispersion of (Jrand = 300kms~^ ( pavis and Peebles 1983| , de Lapparent, 



Geller and Huchra 1988| , Fisher et al. 199^ ). Note that observers know neither Vh nor the true 



distance, but only Vc by observations as described later in this section. 

Absolute magnitudes are assigned to the galaxies according to a galaxian luminosity 
function. Since our maximum likelihood analysis is little affected by the luminosity function, it is 
unnecessary to pursue its true (unknown) shape and amplitude. We adopted a Schechter form 
with = —20.0 mag and a = —1.1, which are representative values found in recent studies 
(Efstathiou, Ellis and Peterson 1988; Loveday et al. 1992; Marzke, Huchra and Geller 1994). 
Then, an intrinsic HI line-profile width W is computed from the absolute magnitude and the 
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TF relation given by equation (^). It is assumed that the TF relation is broadened by its own 
intrinsic scatter with a Gaussian probability function of cjtf = 0.40 mag. Inclination i is assigned 
so that the galaxies have random orientation. Using W, i, and Vc, a projected line-profile width 
W°'"' is obtained via equation (^). Given the intrinsic axial ratio Rq^ = 0.2 , an axial ratio R 
is derived from the inclination using equation @. An apparent magnitude Bj^ is derived from 
the absolute magnitude and the true distance. Internal extinction and K-dimming are computed 
according to a prescription given in RC3. 

Now, these simulated galaxies are observed. The virtual observation imposes observational 
errors to the parameters and to reduce the sample according to a limiting magnitude. First, 
observational errors are added to Vc, W°^'^, R, and Bj^. The errors are assumed to follow 
random Gaussians whose dispersions are estimated from the actual TF sample galaxies. Following 
values are used for the dispersions : AV = lOkms"^ , AW"^'^ = 20kms~^ , Alogi? = 0.06, 
and ABj, = 0.16 mag (§|3.1|). The simulated galaxies are then selected according to a limiting 
magnitude. A selection function S{Bt) is defined as 

r 1 (Bt < 15.0 mag) 

S{Bt) = 1 N{Bt)/No10'^-^^t (15.0 < 5t < 17.0 mag) (16) 
i {Bt> 17.0 mag) , 

where N{B't) is the number of galaxies in the actual TF sample function of Bt (Fig.|l). The 
normalization factor Nq is determined by a least-squares 

fit of iV = TVo • 10°-^^^ to the number 
distribution of the subsample of Bt < 15.0 mag shown in Fig.Q. 

Finally, according to the same selection criteria as imposed on the actual TF sample, i.e., 
Vc > 1000kms-\ i > 45° and 2.30 < log W20 < 2.75 , the simulated magnitude-limited sample is 
reduced to generate a simulated TF sample of 441 galaxies. A velocity-TF distance relation of a 
typical example of the simulated TF samples is shown in Fig.^. 



5.2. Examination of the Performance of the Present Method 

Our TF analysis with the maximum likelihood method is applied to one hundred different 
simulated TF samples. We show in Fig.^ a probability contour map for Hq and utf obtained for 
one of the simulated samples. The contours represent only an internal error associated with the 
likelihood calculation. The input values log/i = —0.15 and ctf = 0.40 mag are reproduced within 
the internal error of the maximum likelihood method (70% confidence level). A distribution of 
the most probable values of {HQ,aTF) is shown in Figjl^ for the one hundred simulated samples. 
Also shown is a distribution of sample averages of log(V^/r7y) for the same samples. In this 
case, the intrinsic scatter aTF is tentatively taken equal to the value obtained by the maximum 
likelihood method. It is clearly demonstrated that the sample average of log{Vc/rTF) is subject 
to a significant systematic error. The .ffo's derived from the maximum likelihood method, on the 
other hand, show a distribution consistent with the input value. 

We further examine the performance of our method by computing a spatial variation of -f^o 
with a simulated TF sample divided into bins of log Vc ■ The intrinsic scatter cttf is fixed for all the 
subsamples in the log Vc bins at 0.41 mag obtained from the analysis shown in Figj9|. The variation 
of Hq is shown in Fig.|ll|a. The evident increase of the sample average of log{Vc/rTF) with log 14 
demonstrates the effect of the sampling incompleteness bias. The apparent increase manifests 
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itself as the systematic error shown in Fig. 10. It is evident, on the other hand, that the maximum 
likehhood method is free from the systematic error caused by the samphng incompleteness bias 
and gives unbiased estimates of Hq's regardless of log Vc- 

As shown in equation (|5|), we corrected for the Malmquist bias assuming an uniform galaxy 
distribution although the actual distribution is inhomogeneous. To examine the effect of ignoring 
the inhomogeneity, we made a comparative simulation using a uniform galaxy distribution. A 
resultant variation of Hq is shown in Fig.^6. Comparison of the two results would reveal the 
effect of inhomogeneity, if any, as a residual. However, we observe no significant difference between 
the two results. Accordingly, we neglect the correction for the inhomogeneous Malmquist bias in 
the present analysis. 



6. RESULT 

6.1. Global Value of Hq 

The maximum likelihood method is applied to our TF sample. The resultant probability 
contour map for Hq and a!^ is shown in Fig.^. Note that the contours represents only the 
internal error of the likelihood calculation. Other error sources are as follows. A zero point 
uncertainty of the TF calibration is evaluated from the standard deviation of the calibrators around 
the TF relation, i.e., 0.45/\/5 = 0.20 mag (§ |4.1| ). The systematic difference of Aq is examined 
among the schemes of BH, RC2, and RSA (Sandage and Tammann 1981), and the difference of Ai 
is examined among RC3, RC2, RSA and Bothun et al. (1985). The result is summarized in Table 
p. Since the same absorption correction schemes are used both for the TF sample and for the local 
calibrators, scheme dependence tends to cancel out. The maximum differences in Aq and in Ai 
are adopted as the respective external errors. These external systematic errors are summarized in 
Table |. 

Finally, we obtain the best estimate of the unbiased global Hq as 



HQ = Q5±2tukms-^Mpc-K (17) 

The first and the second terms of the quoted errors represent the internal random error of the 
maximum likelihood method and the external systematic error, respectively. Note that the latter 
is given by the uppermost and the lowermost boundaries of the linear sum. The errors of Hq are 
often underestimated by ignoring various complicated systematic effects. Our error estimate is 
the result of our best effort to include all the known error sources with conservative error budgets 
assigned to respective sources. The intrinsic scatter of the i?-band TF relation is obtained to be 
a!^ = 0.49 lb 0.05 mag. An analysis with the ten calibrators results in Hq = 68kms~^Mpc~^ 

( Fi) 

with arpp = 0.48 mag. If we use 445 galaxies including four spirals discarded from the TF sample 
(§ [4. 2]) , a!^ slightly increases to 0.54 mag but Hq does not change in practice. 



6.2. Spatial Variation of Hq 

A spatial variation of Hq is examined using the TF sample divided into bins of logV^. 
Dispersion arpp is fixed at 0.49 mag for all the subsamples in the bins; the free parameter is only 
Hq in the maximum likelihood method here. The spatial variation obtained is shown in Fig. 13. 
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It is clearly demonstrated that the sample means of log{Vc/rj,p) are seriously subject to the 
systematic error as Vc increases. On the other hand, the amplitude of the variation of Hq^s is 
consistent with the global Hq within the errors except for poor subsamples at log 14 < 3.6. 



7. ANALYSIS WITH WILLICK'S r-BAND DATA 

7.1. r-band Sample 

The same Hq determination is carried out using r-band CCD data of the Pisces-Perseus 
region given by Willick et al. (1997). The sample (w91pp) is based on UGC and consists of 326 
spirals complete down to niz = 14.5 mag. Application of the same selection criteria as used in our 
TF analysis to w91pp reduced the sample to 271 galaxies. We hereafter refer to this sample of 271 
galaxies as the r-band TF sample. The velocity distribution is shown in Fig.^. The shallower 
limiting magnitude than ours (m^ = 15.5 mag) makes the spatial coverage of this sample less 
extensive than that of the i?-band sample (Fig.^. 



7.2. r-band TF calibration 

The r-band TF relation was calibrated with six local calibrators (Table The same criteria 
as in i?-band leave only two local calibrators (NGC 1365 and 2403). We do not exclude NGC224, 
300, 598 and 3031 simply because the two calibrators are insufficient for reasonable calibration. 
Photometry data for NGC 1365 are taken from ESO-Uppsala catalog ( [Lauberts and Valentijn 



1989) while the data for the others are from Pierce and Tully (1992). Kron-Cousins R and / 
magnitudes are transformed with a color equation r = 0.40 + R — 0.28 (R — I) given by Willick 
(1991). Total magnitudes were corrected for Aq and Ai using BH and Bothun et al. (1985), 
respectively, with a transformation A^^^ = 0.60j4(^). Corrected apparent magnitudes rj, are given 
in Table ^. With the maximum likelihood method, the r-band TF relation is obtained as 

= -20.15 (±0.18) - 7.93 (±1.04) (log W20 - 2.5) . (18) 



The TF diagram of the local calibrators is shown in Fig. 15. The dispersion is 0.42 mag. The slope 
is consistent with that of Pierce & Tully (1992), although we found inconsistency in B (§ [4.11 ). 
This is because the r-band TF diagram in Pierce & Tully (1992) shows no appreciable difference 
in slope between logVF20 > 2.30 and logVF20 < 2.30, which is not the case in the i?-band. Willick 
et al. (1997) gave 7.73 for the slope, which also matchs well with ours. Our correction for Ai 
and for Aq follows the same schemes as theirs, thus the zero point of the TF calibration is also 
comparable with ours. Their zero point is —20.12 {Hq = 65kms~^Mpc~^ assumed), which is in 
good agreement with ours. 

A systematic difference of Alogi?25 = 0.03 is found between our axial ratio R25 and that 



measured by Willick (1991) (Fig. 16). This difference introduces errors of 0.03-0.45 mag in Ai and 
0.3-7 % in W20 depending on galaxy morphology and inclination, respectively. A comparative 
calculation using both sets of the axial ratios for galaxies common between the r- and the i?-band 
samples shows that our axial ratios gives Hq which is 3% smaller than that in r. 



The velocity-distance relation for the r-band TF sample is shown in Fig, 17. We find that 



the r-band sample shows tighter velocity-distance correlation (correlation coefficient = 0.873) 
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than for iiipj (correlation coefficient = 0.804, Fig.^). This may be accounted for by the following 
two reasons. First, the r-band TF relation shows smaller intrinsic scatter than S-band (Pierce 
and Tully 1988). Second, the CCD observations for the r-band sample provide more accurate 
photometry than our photographic photometry. 

Difference between //^ and fii^ is plotted against log Vc in Fig.|l^. There is a systematic 
difference of 0.16 mag. If this is ascribed to the Malmquist bias, i.e., 1.38{a!^ ~^tf ^) ~ ^-^^ mag, 
then we obtain a^p = 0.35 mag from a^p = 0.49 mag (§|6.1]), This value is indeed in agreement 
with what we obtain in the following maximum likelihood analysis (§ |7.3| ). 

Note that the r-band calibrators miss NGC 925 and 4536 which show relatively large deviation 
in B from the calibrated relation (Table §). It is possible that future inclusion of these galaxies 
considerably changes our r-band result. 



7.3. Result 

(r) 

A contour map for Hq and a^p obtained from the maximum likelihood analysis is shown in 
Figjl^. The internal and external errors are assessed in the same way as for the -B-band sample. 
We finally obtain 



i?o = 63 ± 1 km s^^Mpc"^ . (19) 
The quoted errors have the same meaning as equation ( p^ . The intrinsic scatter of the r-band TF 

(r) Tr) 

relation is obtained to be a^p = 0.33 it 0.05 mag. This a>pp well explains the systematic difference 
between ^pp and ^x^pp shown in Fig.p^. 

The spatial variation of Hq is shown in Figj2^. Except for poor subsamples at log Vc < 3.5, 
the values of -f^o's are mostly in agreement with the global value within their errors. Apparent 
increase of a sample mean of log(lc/r^) toward larger log 14 is smaller than that for S-band 
analysis (Fig.^3|). This is consistent with app which is smaller than . 



8. DISCUSSION 

8.1. Global Value of i^o 

We compare our result with those obtained by other studies from the TF relation, the 
period-luminosity relation of Cepheid variables, surface brightness fluctuation (SBF), the planetary 
nebulae luminosity function (PNLF), and supernovae (SNe) type la and E. Extensive reviews 
for these methods are found in Jacoby et al. (1992), van den Bergh (1994), Livio, Donahue and 
Panagia (1997), and Okamura (1998). 



8.1.1. The TuUy-Fisher relation 

IF92 obtained Hq = 78.5 t^!^;? kms-^Mpc-^ using the TF relation calibrated with flve 
calibrators (NGC 224, 598, 300, 2403, 3031). If we adopt the same calibrators as theirs but new 

( Fi) 

Cepheid distances, we obtain a TF calibration of M^p' = -19.44 - 6.19 (log W20 - 2.5). The slope 
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matchs that of IF92, —6.24, but the zero point is 0.28 mag brighter. Using this cahbration, we 
obtain Hq = 74 it 2 ^]^q kms~^Mpc^^. The difference of Ho between IF92 and this calculation is 
fully ascribed to the TF zero point difference of 0.28 mag. 

The Mark III catalog (Willick et al. 1997) provides Hubble velocities of spirals using the TF 
relation. For our galaxies common to the Mark III "w91pp" and "hmcl" samples, we compute 
the Hubble velocities using Hq^^ = 65kms~^Mpc~^ and our TF distance moduli r!^p^ corrected 
for the Malmquist bias. The result of comparison is shown in Fig.^. No systematic difference 
is found in the Hubble velocity, except for a bin 3.5 < logV^B < 3.6 of w91pp. For this bin, we 
observe Hq which is larger than those of other bins (Fig.^). This locally larger value of Hq is 
compatible with the difference of the Hubble velocity in this bin. These facts imply an actual 
infall motion of galaxies in the bin toward the Pisces-Perseus supercluster. 

Recently, Giovanelli et al. (1997a) studied 24 clusters of galaxies out to ~ 9000 km s~^ and 
obtained Hq = 69ib5kms~^Mpc~^. Our result Hq^^ = 65kms~^Mpc~^ is in good agreement with 
theirs. Four of their clusters reside in the Pisces- Perseus region we observed. A comparison of the 
galaxy distances is shown in Fig.|2^. Their distance modulus here is converted from their original 
recession velocity (Giovanelli et al. 1997b) based on Hq = 69 kms~^Mpc~^. The systematic 
difference is —0.06 mag, indicating a satisfactory agreement in the distances. A vertical spread 
is apparent at m — M = 34.4. The spread, A(m — M) ~ 0.5 mag, is consistent with our TF 

/ON 

dispersion = 0.49 ± 0.05 mag. 



8.1.2. Other Methods 

The period-luminosity relation of Cepheid variables, SBF, and PNLF are currently available 
to galaxies up to only a few tens of Mpc. Using Cepheids in galaxies in the Virgo cluster. Pierce 
et al. (1994) and Freedman et al. (1994) obtained = 87 ± 7 and 80 ± 17kms"^Mpc~\ 
respectively. The error quoted by Pierce et al. (1994) did not include the cluster depth effect for 
which Freedman et al. (1994) assumed ±0.35 mag. If the same depth effect is taken into account 
in the result of Pierce et al. (1994), the total error becomes ibl7kms^^Mpc^^. Tanvir et al. 
(1995) identified Cepheid variables in a galaxy in the Leo I group and obtained a smaller value 
-f^o = 69 ± 8kms~^Mpc~^. More recently, Freedman (1998) presented Hq = 73 it 14kms~^Mpc~^ 
in an interim report of the HST key project. From the SBF method with a new zero point 
calibrated with Cepheid observations, Tonry et al. (1997) obtained Hq = 81 it 6kms~^Mpc^^. 
Using the PNLF method, McMillan, Ciardullo and Jacoby (1993) obtained from three galaxies 
in the Fornax cluster Hq = 81 it Skms^^Mpc^^ and 69 it 8kms~^Mpc~^, which depend on the 
choice of the Fornax/Coma distance ratio of 5.25 (Faber et al. 1989) and 6.14 (Aaronson et al. 
1989), respectively. Although PNLF distances are based on the M31 distance of 710 kpc instead 
of 770 kpc which we adopt, the difference reduces their Hq^s by only 3%. 

Using a Hubble diagram of 13 SNe type la calibrated with a single Cepheid distance to 
SN1972E in NGC5253, Riess, Press and Kirshner (1995) obtained Hq = 67 ± 7kms-^Mpc-^ 
within Vc ~ 30000 km s"^ Hamuy et al. (1996) obtained Hq = 63.1 ± 3.4 ± 2.9kms-^Mpc"i in a 
similar spatial extent using a Hubble diagram constructed with 29 SNe of type la and calibrated 
with four SNe (SN 1972E in NGC5253, SN 1937C in IC4182, and SN 1981B in NGC4536, and 
SN 1990N in NGC4639) to which Cepheid distances are available. These Hq^s are both in good 
agreement with ours (see also Saha et al. 1995, Tammann and Sandage 1995, Saha et al. 1996). 
The expanding photosphere method has been developed for the distance determination with type 
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H SNe (Kirshner and Kwan 1974). The method dh'ectly provides the distance to the SNe without 
any external cahbration procedures. Schmidt et al. (1994) obtained Hq = 73 it 6 ib 7kms^^Mpc^^ 
using sixteen SNe type H extending to 14 ~ 15000 km s^^. Their result is slightly larger than but 
in marginal agreement with ours. 

Shanks (1997) compared TF distances to 11 spirals with their Cepheid distances and to 
12 spirals with their SNe la distances. The TF distances were taken from Pierce (1994). They 
concluded that, overall, the TF distance is smaller by 0.43 it 0.11 mag than the Cepheid/SNe 
distances. This leads to that Hq obtained by Pierce (1991) using the Virgo and the Ursa Major 
clusters is revised downwards from Hq = 84 it lOkms^^Mpc^^ to Hq = 69 it 8kms~^Mpc~'^. This 
Hq is in good agreement with ours. 

From the above overview, our Hq = 65 ± 2;^^kms~^Mpc-^ agrees with the results by 
Cepheids observations, the SNe la and H and the Tully-Fisher relation. Although sample galaxies 
are sparse, Cepheids and SNe observations currently provide us with the most reliable distances 
in the local Universe (~ 20 Mpc) and to much further galaxies (~ 500 Mpc), respectively. 
It is a good demonstration that our estimate of the global Hq is consistent with both of these 
results obtained from such different fields of the Universe. It is also convincing that our Hq is 
in agreement with the result of an independent TF study, Giovanelli et al. (1997a), which used 
sparse but widely distributed galaxies. Further comparative studies of the Cepheid, SBF, and 
PNLF distances will work out the current inconsistency of Hq among these three methods. Then, 
Cepheid observations of distant galaxies will give a definite conclusion for the global value of Hq. 



8.2. Bulk Motion 

A bulk motion in the Pisces- Perseus region was investigated by Willick (1990), Han and 
Mould (1992), Courteau et al. (1993), Baffa et al. (1993), and Hudson et al. (1997). Willick 
(1990) claimed that the Pisces-Perseus supercluster and the galaxies behind the supercluster up 
to Vc ~ 8000 km s"^ showed a negative coherent flow. The result was supported by Courteau 
et al.(1993), who argued that all galaxies within a sphere out to at least Vc ~ 6000 km s^^ in 
radius around the Local Group showed parallel streaming flow with amplitude of ~ — 400kms~^ 
toward the Great Attractor. Han and Mould (1992) studied clusters of galaxies lying beyond the 
Pisces-Perseus supercluster, and also detected a negative coherent motion of the clusters. Baffa 
et al. (1993) reported that galaxies around the Pisces-Perseus supercluster are infalling into the 
supercluster, i.e., that the supercluster is shrinking. Most recently, Hudson et al.(1997) studied 
seven clusters in the Pisces-Perseus region and found peculiar motions which are consistent with 
no bulk motion in the CMB frame. They also argued on the basis of x^-statistics that their result 
is also consistent with the negative bulk motion found by Courteau et al. (1993). 



We show in Fig. 23 mean peculiar velocities derived from our Hq and r^y corrected for the 
Malmquist bias. Because of the small number of galaxies at Vc < 3000 km s~^, we discard these 
galaxies from the current discussion. At a farther region, Vc > 8000 km s~^, we find a serious 



effect of the sampling incompleteness bias (see Fig. 13). The same effect is found in Courteau 
et al. (1993) and Mathewson and Ford (1994), but they managed to circumvent the effect by 
comparing the peculiar velocity with that of a control field. As we have no such control data, we 
simply restrict ourselves to the range of 3000 <Vc< 8000 km s^^ for the current discussion on the 
peculiar motion. 

In the i?-band {H^^^ = 65kms^^Mpc^^) we find a mean negative bulk flow of 
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Vp = — 744 lb 92 km s^"*^, while a flow with a smaller amplitude = — 202ib68 kms^^ is obtained 
in r {Hq'^ = 63kms~^Mpc~^). If we change Vc of each galaxies in 3000 < 14 < 8000 km s^^ 
to Vc — Vp^^ and then re-calculate Hq, we obtain Hq = 72kms~^Mpc^^. This value is 
significantly larger than Hq^^ = 65 kms~^Mpc^^. This naive correction of Vc for Vp^^ is probably 
invalid, however, since it substantially reduces the likelihood. The ratio of the likelihood for 
-f^o^^ = 65kms~^Mpc~^ to that for Hq = 72kms~^Mpc~^ is 2.6 x 10~^, which indicates less than 
0.001% probability for the corrected Vc being well fitted with the Hubble flow. A 95% probability 
interval of Vp inferred from the likelihood ratio is —290 < Vp < -|-470kms~ , each of the 
boundary corresponding to Hq = 68 and 61 kms~^Mpc~^, respectively. The inconsistency between 
the 95% probability interval and Vp^^ = — 744kms~^ may be caused by the large uncertainty of 
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the -B-band TF distances due to a large arpp . It also implies that the current i?-band TF analysis 
is not suitable for deducing accurate peculiar motions for individual galaxies. 

In r-band, on the other hand, correction of Vc for Vp = —202 km s leads to 
Hq = 65kms~^Mpc~^ which is not significantly different from Hq~^ = 63 kms~^Mpc~^. Further, 
a 95% probability interval of Vp^^ inferred from the likelihood ratio is —510 < Vp < -|-140kms~^, 

(r) 1 

which is consistent with the mean negative bulk flow Vp = — 202kms^ . Boundaries of the 
95% probability interval of Vp^^ corresponds to Hq = 68 and 62kms~^Mpc~^, respectively. This 
consistent result is convincing for the existence of the negative bulk flow. The range of Vp is 
consistent with those of the earlier results (Table.^). We consider that the r-band TF analysis for 
the peculiar motions seems more reliable than that in S-band. Our result Hq ~ 65kms~^Mpc~^ is 
also consistent with negative bulk flows which were obtained in the previous studies independently 
from an absolute value of Hq. 



9. SUMMARY 

We perform an unbiased estimate of the global Hq using a TF relation with a maximum 
likelihood method. We use a homogeneous, magnitude-limited sample of 441 galaxies 
{niz < 15.5 mag) distributed up to cz ~ 12000 km s^^ in the Pisces-Perseus region. The reliability 
of the method was examined by the Monte-Carlo simulation, and it was shown that the method 
works satisfactorily in deriving an unbiased Hq for the present sample. The results of the present 
analysis are summarized as follows. 

1. We obtain Hq^^ = 65 it 2 kms~^Mpc~^ for the unbiased estimate of the global value of 
the Hubble constant. The first and second error terms represent the internal and external 
errors, respectively. All the known sources of systematic errors are taken into account in 
the estimate of the external error. The same analysis with the r-band TF sample of 271 
galaxies gives -f^o*^^ = 63 it 1 kms~^Mpc~^, which is consistent with H^\ These HQ's 
are in good agreement with IF92 if we allow for the difference of the TF calibration. We 
also find a good agreement for these HQ's with those obtained via Cepheid observations 
for local Universe (Freedman 1998), the Tully-Fisher relation (Giovanelli et al. 1997a) and 
supernovae (Riess, Press and Kirshner 1995, Hamuy et al. 1996) applied to a spatial volume 
larger than ~ 8000 km s~^. 

2. A bulk motion in the Pisces-Perseus region is briefly investigated using our Hq. Analysis of 
the r-band data compiled by Willick et al. (1997) gives a negative motion of 1^ ~ —200 kms~^ 
between the range 3000 < Vc < 8000 km s^^. The amplitude of Vp is consistent with the 
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95% probability interval —510 < Vp < +140 km s^^ inferred from the likelihood ratio. The 
interval is consistent with the bulk flows obtained by Willick (1990), Han and Mould (1992), 
Courteau et al. (1993), and Hudson et al. (1997). Analysis of the S-band data does not give 
a meaningful result because of a larger intrinsic dispersion of the TF relation. 
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., . Number of 

Criteria , 

galaxies 



Photometry sample 




1524 


21cm data available 




957 


Morphological type 


Sa-Sd 


754 


Inclination 


i > 45° 


523 


Recession velocity 


Vc > lOOOkms^i 


514 


21cm line-profile width 


2.30 < log W20 < 2.75 


W449 


Absolute magnitude residual 


{M^p' -Mb\< 1.62 mag 


(fa) 441 



Table 1: Sample selection for the Tully-Fisher analysis 



Note. — (a) Spiral sample, (b) TF sample. 



Name 
(1) 


Type 
(2) 


M 

[mag] 

(3) 




[mag] 
(4) 


[mag 
(5) 


[kms-i" 
(6) 


[kms~-^ 
(7) 


logi?25 
(8) 


Mb - M^p^ 

[mag] 
(9) (10) 


N224t (M31) 


.SAS3.. 


(a) 24.44 ( 


13) 


3.34 (.02) 


2.67 


536 


-300 


0.49 


0.13 


-0.01 


N598t (M33) 


.SAS6.. 


(^)24.64 ( 


09) 


5.74 (.03) 


5.49 


199 


-179 


0.23 


0.11 


0.09 


N300 


.SAS 7.. 


('=)26.66 ( 


10) 


8.53 (.05) 


7.89 


150 


142 


0.15 


0.48 


0.49 


N2403 


.SXS6.. 


('^)27.51 ( 


24) 


8.41 (.07) 


7.91 


244 


131 


0.25 


0.39 


0.34 


N30311' (M81) 


.SAS 2.. 


(^)27.80 ( 


20) 


7.40 (.03) 


6.36 


434 


-34 


0.28 


0.55 


0.42 


N925 


.SXS 7.. 


(•^) 29.84 ( 


16) 


10.00 (.11) 




215 


553 


0.25 


-0.69 


-0.72 


N3368 (M96) 


.SXT2.. 


(9)30.32 ( 


16) 


9.86 (.13) 




354 


897 


0.16 


0.38 


0.26 


N4536 


.SXT4.. 


W3I.IO ( 


13) 


10.56 (.08) 




336 


1804 


0.37 


-0.45 


-0.53 


N1365t 


.SBS3.. 


W3I.32 ( 


17) 


9.90 (.07) 


8.94 


398 


1662 


0.26 


-0.67 


-0.79 


N4639t 


.SXT4.. 


(^■)32.00 ( 


23) 


11.92 (.10) 




356 


1010 


0.17 


0.72 


0.60 



Table 2: Local calibrators 



Note. — Col(2): Morphological type taken from RC3. Col(3): Distance modulus and its error, /i = 18.50 mag is 
assumed for LMC. References are: (a) Freedman and Madore (1990), (6) Freedman, Wilson and Madore (1991), (c) 
Freedman et al. (1992), (d) Freedman and Madore (1988), (e) Freedman et al. (1994), (/) Silbermann et al. (1996), 
{g) Tanvir et al. (1995), (h) Saha et al. (1996), (i) Freedman, Madore and Kennicutt (1997), and {j) Sandagc ct 
al. (1997). Col(4): _B-band magnitude corrected for the Galactic and the internal extinction and its error taken 
from RC3. Col(5): r-band magnitude corrected for the Galactic and the internal extinction prescribed by Willick 
et al. (1997). Col(6): Observed 21cm line-profile width taken from RC3. Col(7): Heliocentric recession velocity 
taken from RC3. Col(8): Logarithmic major-to- minor axial ratio taken from RC3. Col(9),(10): Difference between 
the B-band absolute magnitude and that derived from the Tully-Fisher relation calibrated with five (col.9) and ten 
(col.lO) calibrators. ^ See text. 
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Ag 






A, 




RC2-BH 


RSA-BH 


RC2-RC3 


RSA-RC3 Botliun-RC3 


TF sample 


0.14 (.06) 


-0.07 (.07) 


-0.12 (.07) 


0.36 (.15) -0.17 (.06) 


Local calibrators 


0.15 (.05) 


-0.06 (.03) 


-0.14 (.07) 


0.31 (.14) -0.15 (.03) 


Difference 


-0.01 


-0.01 


0.02 


0.05 -0.02 



Dispersion is given in parentheses. 



Table 3: Difference of average extinction corrections for the TF sample 



o Internal random error 
Statistics (Fig.|l|) (±2kms-iMpc-i) ±2.7% 

o External systematic error 



Tully-Fisher calibration (±0.20 mag) ±9.7% 

LMC distance (±0.1 mag) ±4.6% 

Photometric zero point (±0.05 mag) ±2.3% 

Axial ratio (A log i?25 = 0.03) ±9% 

Galactic extinction (±0.01 mag) ±0.5% 

Internal extinction (±0.05 mag) ±2.3% 

Parameters {(Jrand, ^W^^' , Alogi;25, A^r) (Alog/t = 0.005) ±1.2% 

(Arithmetic sum) ±30%, -21% 



Table 4: Summary of the errors for Hq 



Reference 




R.A. 


Volunie 
(kms~ ) 


(kms ) 


This work (B, i?o = 68 - 61 km s 


-iMpc-i) 


22^*- 




3000-8000 


(-290, ±470) 


(r, Hq = eSkms^^Mp 




22^*- 




3000-8000 


-202 ± 68 


(r, iJo = 68-62kms" 


^Mpc-i) 


22'*- 




3000-8000 


(-510, ±140) 


WiUick (1990) 




22'*- 




3800-6000 


-441 ± 49 


Han and Mould (1992) 




o'*- 


-2'* 


4000-7500 


-420 ± 115 


Courteau et al. (1993) 




22'*- 




4000-9000 


-383 ± 58 


Hudson et al. (1997) 




o'*- 


-4'* 


4000-6000 


-60 ± 220 



Table 5: Comparison with previous works of the bulk flow in the Pisces-Perseus region 



- 21 - 



REFERENCES 

Aaronson, M., Bothun, G., Mould, J., Huchra, J., Schommer, R. A., Cornell, M. E. 1986, ApJ, 
302, 536 

Aaronson, M., Bothun, G., Cornell, M. E., Dawe, J. A., Dickens, R. J., Hall, P. J., Han, M., 
Huchra, J., Lucey, J. R., Mould, J. R., Murray, J. D., Schommer, R. A., Wright, A. E. 
1989, ApJ, 338, 654 

Baffa, C, Chincarini, G., Henry, R. B. C, Manousoyanaki, J. 1993, A&A, 280, 20 

Bernstein, G. M., Guhathakurta, P., Raychaudhury, S., Giovanelli, R., Haynes, M. P., Herter, T., 
Vogt, N. P. 1994, AJ, 107, 1962 

Bothun, G. D., Aaronson, M., Schommer, B., Mould, J. R., Huchra, J., Sullivan HI, W. T. 1985, 
ApJS, 57, 423 

Bottinelli, L., Gouguenheim, L., Paturel, G., de Vaucouleurs, G. 1983, A&A, 118, 4 

Bottinelli, L., Fouquc, P., Gouguenheim, L., Paturel, G., Teerikorpi, P. 1987, A&A, 181, 1 

Bureau, M., Mould, J. R., Staveley-Smith, L. 1996, ApJ, 463, 60 

Cornell, M. E., Aaronson, M., Bothun, G. D., Mould, J. R. 1987, ApJS, 64, 507 

Courteau, S., Faber, S. M., Dressier, A., Willick, J. A. 1993, ApJ, 412, L51 

Davis, M., Peebles, P. J. E. 1983, ApJ, 267, 465 

de Lapparent, V., Geller, M. J., Huchra, J. P. 1988, ApJ, 332, 44 

de Vaucouleurs, G., dc Vaucouleurs, A., Corwin, H. G., Jr. 1976, Second Reference Catalogue of 
Bright Galaxies (RC2) , (Austin: University of Texas Press) 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr, Bura, R. J., Paturel, G., Fouque, P. 
1991, Third Reference Catalogue of Bright Galaxies (RC3), (New York: Splinger-Verlag) 

Faber, S. M., Wegner, G., Burstein, D., Davies, R. L., Dressier, A., Lynden-Bell, D., Terlevich, R. 
J. 1989, ApJS, 69, 763 

Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., Huchra, J. P. 1994, MNRAS, 267, 927 

Fouque, P., Bottinelli, L., Gouguenheim, L., Paturel, G. 1990, ApJ, 349, 1 

Freedman, W. L. 1998, in Proc. lAU symp. 183 on Cosmological Parameters and Evolution of the 

Universe, in press 

Freedman, W. L., Madore, B. F. 1988, ApJ, 332, L63 
Freedman, W. L., Madore, B. F. 1990, ApJ, 365, 186 

Freedman, W. L., Madore, B. F., Kennicutt, R. C. 1997, in The Extragalactic Distance Scale 

(eds. Livio, M., Donahue, M., Panagia, N.), STScI Symp. Ser. 10 (Cambridge: Cambridge 
University Press) 

Freedman, W. L., Wilson, C. D., Madore, B. F. 1991, ApJ, 372, 455 



-22- 



Freedman, W. L., Madore, B. F., Hawley, S. L., Horowitz, I. K., Mould, J., Navarrete, M., 
Sallmen, S. 1992, ApJ, 396, 80 

Freedman, ,W. L., Madore, B. F., Mould, J. R., Ferrarese, L., Hill, R., Kennicutt, R. C, Saha, A., 
Stetson, P., Graham, J. A., Ford, H., Hoessel, J. G., Huchra, J., Hughes, S. M., Illingworth, 
G. D. 1994, Nature, 371, 757 

Freudling, W., Da Costa, L. N., Wegner, G., Giovanelh, R., Haynes, M. P., Salzer, J. J. 1995, AJ, 
110, 920 

Freudling, W., Martel, H., Haynes, M. P. 1991, ApJ, 377, 349 

Fukugita, M., Okamura, S., Tarusawa, K., Rood, H.J., Williams, B. A. 1991, ApJ, 376, 8 
Fukugita, M., Okamura, S., Yasuda, N. 1993, ApJ, 412, L13 
Giovanelli, R., Haynes, M. P. 1985, AJ, 90, 2445 
Giovanelli, R., Haynes, M. P. 1989, AJ, 97, 633 

Giovanelli, R., Haynes, M. P., Myers, S. T., Roth, J. 1986, AJ, 92, 250 

Giovanelli, R., Haynes, M. P., Salzer, J. J., Wegner, G., Da Costa, L. N., Freudling, W. 1994, AJ, 
107, 2036 

Giovanelli, R., Haynes, M. P., Da Costa, L. N., Freudling, W., Salzer, J. J., Wegner, G. 1997a, 
ApJ, 477, LI 

Giovanelli, R., Haynes, M. P., Herter, T., Vogt, N. P., Wegner, G., Salzer, J. J., Da Costa, L. N., 
Freudling, W. 1997b, AJ, 113, 22 

Hamabe, M., Ichikawa, T. 1991, Astronomical Data Analysis Software and Systems I, ASP Conf. 
Ser., vol.25 (eds. Worrall et al.) (San Franscisco: ASP), 325 

Hamuy, M., Phillips, M. M., Maza, J., Suntzeff, N. B., Schommer, R. A., Aviles, R. 1995, AJ, 109, 
1 

Han, M. 1992, ApJS, 81, 35 
Han, M., Mould, J. R. 1990, ApJ, 360, 448 
Han, M., Mould, J. R. 1992, ApJ, 396, 453 
Haynes, M. P., Giovanelh, R. 1986, ApJ, 306, L55 



Hudson, M. J., Lucey, J. R., Smith, R. J., Steel, J. astro-ph /9707026 



Ichikawa, S., Okamura, S., Watanabe, M., Hamabe, M., Aoki, T., Kodaira, K. 1987, Ann. Tokyo. 
Astron. Obs. 2nd Ser., 21, 285 

Ichikawa, T., Fukugita, M. 1992, ApJ, 394, 61 (IF92) 

Jacoby, G. H., Branch, D., Ciardullo, R., Davies, R. L., Harris, W. E., Pierce, M. J., Pritchet, C. 
J., Tonry, J. L., Welch, D. L. 1992, PASP, 104, 599 

Kirshner, R. P., Kwan, J. 1974, ApJ, 193, 27 



-22,- 



Kodaira, K., Okamura, S., Ichikawa, S., Hamabe, M., Watanabe, M. 1990, Photometric Atlas of 
Northern Bright Galaxies, (Tokyo: University of Tokyo Press) 

Kolatt, T., Dckcl, A., Ganon, G., Willick, J. A. 1996, ApJ, 458, 419 

Landolt, A. U. 1992, AJ, 104, 340 

Lauberts, A., Valentijn, E. A. 1989, The Surface Photometry Catalogue of the ESO-Uppsala 
Galaxies, (Munich: European Southern Observatory) 

Livio, M., Donahue, M., Panagia, N., eds 1997, The Extragalactic Distance Scale, STScI Symp. 
Ser. 10 (Cambridge: Cambridge University Press) 

Longo, G., de Vaucouleurs, A. 1983, A General Catalogue of Photoelectric Magnitudes and Colors 
in the U,B,V System of 3,578 Galaxies Brighter Than the 16-th V-magnitude (1936-1982), 
(Austin: University of Texas) 

Loveday, J., Peterson, B. A., Efstathiou, G., Maddox, S. J. 1992, ApJ, 390, 338 

Madore, B. F., Preedman, W. L. 1991, PASP, 103, 933 

Malmquist, K. G. 1922, Med.Lunds.Astr.Obs, 16, No.23 

Marzke, R. O., Huchra, J. P., Geller, M. J. 1994, ApJ, 428, 43 

Mathewson, D. S., Ford, V. L. 1994, ApJ, 434, L39 

McMillan, R., Ciardullo, R., Jacoby, G. H. 1993, ApJ, 416, 62 

Nilson, P. 1973, Uppsala General Catalogue of Galaxies (UGC), (Uppsala: Uppsala Obs. Ann.) 

Okamura, S. 1998, in Proc. The 11th Nishinomiya- Yukawa Memorial Symposium on 'Physics in 
the 21st Century' 

Pierce, M. J. 1991, in Observational Tests of Cosmological Inflation (eds. Shanks, T., Banday, A., 
Ellis, R. S., Frenk, C. S., Wolfendale, A. W.) (Dordrecht: Kluwer) 

Pierce, M. J. 1994, ApJ, 430, 53 

Pierce, M. J., Tully, R. B. 1988, ApJ, 330, 579 

Pierce, M. J., Tully, R. B. 1992, ApJ, 387, 47 

Pierce, M. J., Welch, D. L., McClure, R. D., van den Bergh, S., Racine, R., Stetson, P. B. 1994, 
Nature, 371, 385 

Riess, A. G., Press, W. H., Kirshner, R. P. 1995, ApJ, 438, L17 

Saha, A., Sandage, A., Labhardt, L., Schwengeler, H., Tammann, G. A., Panagia, N., Macchetto, 
F. D. 1995, ApJ, 438, 8 

Saha, A., Sandage, A., Labhardt, L., Tammann, G. A., Macchetto, F. D., Panagia, N. 1996, ApJ, 
466, 55 

Sandage, A., Tammann, G. A. 1981, A Revised Shapley-Ames Catalog of Bright Galaxies (RSA), 
(Washington: Carnegie Inst, of Washington Pub.) 



- 24 - 



Sandage, A. et al. 1997, ApJ submitted. 
Schechter,P.L., 1980, AJ, 85, 801 

Schmidt, B. P., Kirshner, R. P., Eastman, R. G., Phillips, M. M, Suntzeff, N. B., Hamuy, M., 
Maza, J., Aviles, R. 1994, ApJ, 432, 42 

Shanks, T. 1997, MNRAS, 290, L77 

Silbermann, N. A., Harding, P., Madore, B. F., Kennicutt, R. C, Saha, A., Stetson, P. B., 
Preedman, W. L., Mould, J. R., Graham, J. A., Hill, R. J., Turner, A., Bresohn, F., 
Ferrarese, L., Ford, H., Hocsscl, J. G., Han, M., Huchra, J., Hughes, S. M., Illingworth, G. 
D., Phelps, R., Sakai, S. 1996, ApJ, 470, 1 

Smoot, G. F., Bennett, C. L., Kogut, A., Aymon, J., Backus, C., de Amici, G., Galuk, K., 

Jackson, P. D., Kcegstra, P., Rokke, L., Tenorio, L., Torres, S., Gulkis, S., Hauser, M. G., 
Janssen, M. A., Mather, J. C., Weiss, R., Wilkinson, D. T., Wright, E. L., Boggess, N. W., 
Cheng, E. S., Kelsah, T., Lubin, P., Meyer, S., Moseley, S. H., Murdock, T. L., Shafer, R. 
A., Silverberg, R. F. 1991, ApJ, 371, LI 

Strauss, M. A., Willick, J. A. 1995, Phys. Rep., 261, 271 

Tammann, G. A., Sandage, A. 1995, ApJ, 452, 16 

Tanvir, N. R., Shanks, T., Ferguson, H. C., Robinson, D. R. T. 1995, Nature, 377, 27 

Tonry, J. L., Blakeslee, J. P., Ajhar, E. A., Dressier, A. 1997, ApJ, 475, 399 

Tully, R. B., Fisher, J. R. 1977, A&A, 54, 661 

van den Bergh, S. 1994, PASP, 106, 1113 

Watanabe, M. 1996, Ph.D. thesis. University of Tokyo 

Willick, J. A. 1990, ApJ, 351, L5 

Willick, J. A. 1991, Ph.D. thesis. University of California 

Wilhck, J. A., Courteau, S., Faber, S. M., Burstein, D., Dekel, A., Strauss, M. A. 1997, ApJS, 
109, 333 

Yasuda, N., Okamura, S., Fukugita, M. 1995, ApJS, 96, 359 

Zwicky, F., Herzog, E., Wild, P., Karpowicz, M., Kowal, C. T. 1961-68, Catalogue of Galaxies and 
of Clusters of Galaxies (CGCG), (Pasadena: California Institute of Technology) 



This preprint was prepared with the AAS macros v4.0. 



- 25 - 




1.0 - 



O 

o 



o 



0.5 - 



0.0 



1 1 1 1 1 1 1 1 1 1 1 
(b) 


1 'y' 






F}4 








□ □ ° 












,'' 1 1 1 1 1 1 1 1 1 1 1 


1 . 



0.0 



0.5 

logR25(RC3) 



1.0 



Fig. 1. — Comparisons (a) of the total magnitude Bt and (6) of the logarithmic isophotal major- 
to-minor axial ratio logi225j between RC3 data and ours. Filled circles represent galaxies for which 
RC3 data are based on surface photometry. A encircled cross indicates PGC5473 (see text). 
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Fig. 2. — Comparisons (a) of the isophotal magnitude B25 and (6) of the logarithmic isophotal 
major-to-minor axial ratio logi?25, between the data of Cornell et al (1987) and ours. Encircled 
crosses indicate PGC5473 and PGC3133 (see text). 
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Fig. 3. — Velocity distribution of the spiral sample (A'^ = 449). The velocity is expressed in the 
CMB-rest frame. 




Fig. 4. — A histogram of the total magnitude Bt for the spiral sample {N = 449, thick line) and 
the photometry sample {N = 1524, thin line). A dashed line shows a number ratio of the former 
sample to the latter. 
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Fig. 5. — The i?-band Tully-Fisher diagram for the local calibrators given in Table 2. The NGC 
number is labeled. Open circles represent five galaxies discarded from the fiducial calibration (see 
text). The solid and dashed lines show regression lines for the fiducial set of five calibrators and 
for all the ten calibrators, respectively. 
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Fig. 6. — Difference between tlie TF magnitude Alj^p and tlie magnitude from the Hubble's 
law witli log/i = —0.15, plotted against logi?25 for the spiral sample (N = 449). Crosses represent 
galaxies whose images are contaminated by other images. Upper and lower dashed lines delineate 
the boundaries \Mj.p' — Mb\ = 1.62 mag. Eight galaxies (encircled crosses and triangles) are out of 
the boundaries and hence discarded from the TF analysis. The enclosed crosses indicate galaxies 
with unreliable observational data (see text). 
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Fig. 7. — The velocity-TF distance relation for the spiral sample {N = 449) . Symbols are the same 
as in Fig.6. A dashed line represents an arbitrary Hubble's law with log/t = —0.15. 
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Fig. 8. — Same as Fig.7 but for one of the simulated TF samples {N = 441). 
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Fig. 9. — Probability contour map for Hq and aTF obtained from one of the simulated TF samples. 
Confidence levels at 70%, 95% and 99% are shown. Horizontal and vertical dashed lines indicate 
the input values log/i = —0.15 and axF = 0.40 mag, respectively, which are assumed a priori for 
the simulated sample. 
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Fig. 10. — The most probable values of Hq and cttf for one hundred different simulated TF 
samples. Filled squares represent those obtained by the maximum likelihood method, while open 
circles represent sample averages of log(V^/rTF) for log h with gtf a priori set equal to that obtained 
with the maximum likelihood method. 
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Fig. 11. — Spatial variation of Hq obtained (a) from the simulated TF sample shown in Fig.9 and 
(6) from a simulated TF sample with an uniform spatial distribution. The intrinsic scatter axF 
is fixed at 0.41 mag for all the subsamples in logV^ bins. Dots represent the individual galaxies, 
to which the Hubble ratio Vc/ttf is given as Hq. Filled squares represent the most probable Hq 
in each bin obtained by the maximum likelihood method, while open circles represent the sample 
average of log{Vc/rTF)- The squares and the circles are shifted from the center of each log Vc bin to 
avoid the overlap. Error bars are taken from the 70% error obtained for each subsample. A dashed 
line indicates the input value log/i = — 0.15. 
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Fig. 12. — Probability contour map for Hq and arpp for the TF sample of 441 galaxies. A plus 
indicates the most probable values. Contours are drawn at 70%, 95% and 99% confidence levels. A 
triangle and a square indicate the results obtained from the ten calibrators and from the inclusion 
of four deviant galaxies, respectively (see text). 
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Fig. 13. — Spatial variation of Hq obtained from the TF sample. The intrinsic scatter is fixed 
at 0.49 mag for all the subsamples in log Vc bins. Dots represent the individual galaxies, to which 
the Hubble ratio Vc/r^ is given as Hq. Filled squares represent the most probable Hq in each bin 
obtained by the maximum likelihood method, while open circles represent the sample average of 
log{Vc/rj,p ). The squares and the circles are shifted from the center of each log Vc bin to avoid the 
overlap. Error bars are taken from the 70% error obtained for each subsample. Two dashed lines 
indicate the 70% confidence levels Hq = 63 and 67kms~^Mpc~^ shown in Fig. 12. 
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Fig. 14. — Velocity distribution of the r-band TF sample {N = 271). The velocity is expressed in 
the CMB-rest frame. 




Fig. 15. — The r-band Tully-Fisher diagram for the local calibrators given in Table 2. The NGC 
number is labeled. The regression line shows equation (18). 
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Fig. 16. — A comparison of the major-to-minor axial ratios between the data of WiUick (1991) and 
ours. 
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Fig. 17. — The velocity-TF distance relation for the r-band TF sample. A dashed line represents 
an arbitrary Bubble's law with log/i = —0.15. 
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Fig. 18. — Difference between the TF distance moduli Hrpp and Hrpp. 




Fig. 19. — A probability contour map for Hq and aj,p obtained from the r-band TF sample. 
Confidence levels at 70%, 95% and 99% are shown. 
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Fig. 20. — Spatial variation of Hq obtained from the r-band TF sample. The intrinsic scatter 

(r) 

arpp is fixed for all the subsamples in logV^ bins to 0.33 mag obtained in Fig. 19. Dots represent 

, (r) . 

the individual galaxies, to which the Hubble ratio Vc/rj,p is given as Hq. Filled squares represent 
the most probable Hq in each bin obtained by the maximum likelihood method, while open circles 
represent the sample average of log(K/r^). The squares and the circles are shifted from the center 
of each log Vc bin to avoid the overlap. Error bars are taken from the 70% error obtained for each 
subsample. Two dashed lines indicate the 70% confidence levels Hq = 62 and 64kms~^Mpc~^ 
shown in Fig. 19. 
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Fig. 21. — Relative difference of the Hubble velocity between Mark IE catalog (Willick et al. 1997) 
and ours, plotted against the Mark HI Hubble velocity, (a) "w91pp" (b) "hmcl". Dots represent 
individual galaxies and solid squares show median in each log Vm3 bin. Error bars are calculated 
as a standard deviation divided by a square root of the number of galaxies in each bin. 
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Fig. 22. — Comparison of the distance modulus between Giovanelli et al. (1997b) and this study. 
Squares, crosses, closed and open triangles represent galaxies in NGC 383 group, NGC 507 group, 
A 262, and A2634/A2666, respectively. Distance moduli for their galaxies are based on Hq = 69 
kms^^Mpc""^ from Giovanelli et al. (1997a). 
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Fig. 23. — Peculiar velocities plotted against the CMB-rest recession velocity in B (upper) and r 
(lower) . A solid horizontal line shows a mean value {Vi'^^ = 744kms-i in B and vi""^ = 202 km s'^ 
in r) in the range 3000 < K < 8000 km 8"^ 



